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We construct exact ground states of interacting electrons on triangle and diamond Hub- 
bard chains. The construction requires (i) a rewriting of the Hamiltonian into positive 
semidefinite form, (ii) the construction of a many-electron ground state of this Hamiltonian, 
and (iii) the proof of the uniqueness of the ground state. This approach works in any dimen- 
sion, requires no integrability of the model, and only demands sufficiently many microscopic 
parameters in the Hamiltonian which have to fulfill certain relations. The scheme is first 
employed to construct exact ground state for the diamond Hubbard chain in a magnetic 
field. These ground states are found to exhibit a wide range of properties such as flat-band 
ferromagnetism and correlation induced metallic, half-metallic or insulating behavior, which 
can be tuned by changing the magnetic flux, local potentials, or electron density. Detailed 
proofs of the uniqueness of the ground states are presented. By the same technique exact 
ground states are constructed for triangle Hubbard chains and a one-dimensional periodic 
Anderson model with nearest-neighbor hybridization. They permit direct comparison with 
results obtained by variational techniques for /-electron ferromagnetism due to a flat band 
in CeRh 3 B 2 . 

§1. Introduction 

Electronic correlations in solids can result in a high sensitivity of these systems to 
small changes in external parameters such as temperature, pressure, magnetic field, 
or doping. This may then lead to strongly non- linear responses such as large 
resistivity and volume changes, a strong thermoelectric response, colossal magnetore- 
sistance, and high-T c superconductivity. The complexity displayed by electronically 
correlated materials is not only the topic of fascinating fundamental research, but 
is also significant for technological applications, e.g., in sensors, switches, and in 
spintronics. 

Even highly simplified models of interacting electrons are known to display a 
dazzling diversity of phases, e.g., metallic, insulating, and many types of long-range 
ordered states. 4 )' 5 ) Theoretical investigations of correlation effects almost always 
make some kind of approximation which cannot, in general, explain the full com- 
plexity of a model in detail. In this situation exact solutions of correlation models 
are particularly important and useful. 

In this paper we discuss a general strategy for the construction of exact many- 
electron ground states of lattice models such as the Hubbard model. The approach 
is based on three distinct steps: (i) the rewriting of the Hamiltonian in positive 
semidefinite form, (ii) the construction of exact ground states for this Hamiltonian, 
and (iii) the proof of the uniqueness of the ground states constructed thereby. 
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This strategy was recently employed to construct a class of exact ground states 
of a three-dimensional periodic Anderson model (PAM), including the conventional 
PAM on regular Bravais lattices at and above 3/4 filling, 6 )' 7 ) and to the diamond 
Hubbard chain in a magnetic field. 8 ) In both wide range of features such as 

metallic, non-Fermi liquid, and insulating phases in the case of the PAM, 6 )' 7 ) and 
flat-band ferromagnetism, correlation-induced metallic, half-metallic, or insulating 
behavior in the case of the diamond Hubbard chain 8 ) was found. These exact results 
prove that the phase diagram of even seemingly simple quantum mechanical models 
of many-particle systems can indeed be surprisingly complex. 

The construction of exact ground states by casting the Hamiltonian into pos- 
itive semidefinite form was originally applied by Brand and Giesekus to the Hub- 
bard model and the PAM on special perovskite-type decorated lattices in dimen- 
sions d > 2 at U = oo. 9 ) The procedure was generalized by Strack 10 ) to the PAM 
and the extended Emery model on regular lattices in low dimensions, 11 ) and to the 
extended Hubbard model. The method was subsequently used to derive rigorous 
criteria for the stability of saturated ferromagnetism in general Hubbard-type mod- 
els. 12 ) The possibility of superconductivity in the extended Hubbard model and its 
strong coupling limit was also investigated. 13 )~ 15 ) Furthermore, next-nearest neigh- 
bor terms in the Hamiltonian were shown to cause charge density waves, degenerate 
non-magnetic localized phases, and phase separation. 16 ) Introducing operators to 
rewrite the Hamiltonian in positive semidefinite form which are defined on the entire 
unit cell of a d-dimensional lattice (e.g., on plaquettes in d = 2 or cubes in d = 3) 
rather than only on bonds, Orlik and Gulacsi 17 ) investigated the PAM on regu- 
lar two- and three-dimensional lattices at U = oo, finding insulating non-magnetic 
phases. This method was used by Sarasua and Continentino 18 ) to analyze the one- 
dimensional PAM in more detail. Gurin and Gulacsi 19 ) and Gulacsi 20 ) demonstrated 
that it is also possible to investigate the PAM in d = 2, 3 at U < oo . The exact 
ground states constructed by Gulacsi and Vollhardt for the PAM at 3/4 6 ) and 1/4 7 ) 
gave the first mathematical proof for the presence of ferromagnetism in this model in 
dimension d = 3. These authors also detected a strong anomaly in the compressibil- 
ity of the system at the boundary of the parameter region where the non-magnetic 
localized phase is stable — a feature which is typical for / electron systems such 
as Cerium. In this connection the uniqueness of the constructed states was proved 
explicitly for the first time. 

Most recently, the present authors constructed exact ground states of interacting 
electrons on the diamond Hubbard chain in a magnetic field. 8 ) This model was 
shown to exhibit a wide range of properties such as flat-band ferromagnetism and 
correlation induced metallic, half-metallic or insulating behavior which can be tuned 
by changing the magnetic flux, local potentials, or the electron density. 

In this paper we present details of the construction of exact many-electron 
ground states on the diamond Hubbard chain and of their physical properties. 8 ) 
In particular, in the Appendices A and B we provide details of the mathemati- 
cal proof of the uniqueness of the ground-state solutions for this model. We also 
construct exact ground states on triangle Hubard chains and show that within a 
particular subspace of parameters the triangle Hubbard chain may be mapped to 
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a one-dimensional PAM. The exact solutions obtained for this model can be ap- 
plied to investigate the correlated electron material CeRh3E>2, for which Kono and 
Kuramoto 21 ) had constructed an approximate solution. 

The paper is structured as follows: In Sec. 2 we discuss the general strategy 
for the construction of exact many-electron ground states. Exact ground states are 
constructed for the diamond Hubbard chain in Sec. 3. Details of the proof of the 
uniqueness of these ground states are presented in the appendices A and B. The 
technique is employed in Sec. 4 to obtain exact ground states also for the triangle 
Hubbard chain, and the results are applied to explain properties of CeRli3B2. 

§2. Construction of exact many-electron ground states: General 

Strategy 

Our construction of exact many-electron ground states is based on the following 
three steps. 

2.1. Step 1: Rewriting the Hamiltonian in positive semidefinite form 

In the first step the many-electron Hamiltonian is cast into positive semidefinite 
form. This means that one rewrites the Hamiltonian in terms of a number of positive 
semidefinite operators P n as 

L 

H = H + H u = Y,Pn + E g = H' + E g , (2-1) 

n=l 

where E g is the ground-state energy. Here (2-1) is supposed to be an operator 
identity, involving no approximations. The positive semidefinite nature of P n is 
expressed by the relation (&\P n \ty) > 0, which must be satisfied for all Due to 
this inequality, the minimal eigenvalue of P n is zero. This property, i.e., the existence 
of a well-defined lower bound, namely zero, in the spectrum of H' = H — E g in (2-1) 
is the corner stone of the procedure. Simple examples for a positive semidefinite 
operator are P n = &f2, or P n = or P = £\ P\, where the positive semidefinite 

operator P\ = n-^h^ — (n^j + n^jj + 1 assumes its minimum eigenvalue if at least 
one electron is present on the site i. 

The scalar E g in (2-1), i.e., the ground state energy, is usually a complicated 
function of the parameters entering in H. It should be noted that the decomposition 
of H in positive semidefinite form is, in general, not unique and can be obtained in 
different ways. Each transformation is valid only in certain parts 1Z P of the parameter 
space. 

2.2. Step 2: Construction of the exact many- electron ground states 

After the transformation of a given Hamiltonian H into a particular positive 
semidefinite form valid in a parameter space region 1Z P , the exact ground state of 
H needs to be constructed. Since H' = J2n=i^n, and the minimal eigenvalue of 
each P n is zero, this requires the construction of the most general state \\P g ) which 
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satisfies 

PnWg) = 0, (2-2) 

for all n. Consequently, 

H\9 g ) = E g \V g ) (2-3) 

holds, where \& g ) is the ground state, with E g as the ground state energy in 1Z P . The 
actual construction of \SP g ) depends on the explicit expression of the operators P n in 
(2-1). Therefore each case must be individually analyzed. Nevertheless the following 
guideline may be applied (here we restrict ourselves to fermionic particles): 

i) For operators P n = Wfi the ground state has, in general, the form \<fi g ) = 
Ot|0), where |0) is the vacuum, and the anti-commutation relation O'} = 
holds. Indeed, in this case one has P n \<j>g) = 0. To construct the global ground state 
one then tries to extend this property to all operators P n . 

ii) Likewise, for operators P n = QW the ground state construction makes use 
of (12^) 2 = and starts from states of the form \tp g ) = /F|0), s ince in this case 
Pntyg) = 0- To obtain \^ g ) from \ip g ) the property P n \ip g ) = 0, valid for a particular 
n, must be extended to all operators P n . 

2.3. Step 3: Proof of the uniqueness of the ground states 

Finally the constructed ground state has to be shown to be unique, i.e., one 
must prove that there exist no other linearly independent eigenstates of H in 1Z P 
with the same energy E g . This step employs the property of the kernel of H' which 
is defined as the Hilbert subspace of states |^} with H'\<P) = 0, e.g. 

ker(H') := { |<f>) | H'\$) = 0}. (2-4) 

Since H' = Yln=iPn, one nas ker(H') = C\n=i^ er (Pn)- To prove the uniqueness 
one must show that \\P g ) spans ker(H'). For this it is necessary and sufficient to 
prove that 

(i) \& g ) is an element of the intersection of the kernels of all positive semidefinite 
operators P n , i.e., \\P g ) G ker(H') = Dn=i ker(P n ), and that 

(ii) all states |#) G f|n=i ker(P n ) can be written in terms of the constructed ground 
state \& g ). 

We note that if condition (ii) is not satisfied the constructed ground state is not 
unique, implying that other, linearly independent states with the same lowest energy 
E g exist. The uniqueness of the ground states constructed in this paper is proved in 
Appendices A and B. 

2.4. Discussion 

The procedure described above works, in principle, in any dimension and for 
any model, independent of the concept of integrability. On the other hand, a decom- 
position of H into positive semidefinite form and a subsequent construction of the 
ground state may not always be possible — especially for Hamiltonians with a very 
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simple structure. Indeed, the procedure works the better the more microscopic pa- 
rameters (e.g., hopping amplitudes, bands, etc.) enter in H. This is due to the fact 
that in a high-dimensional parameter space it is easier to find a subspace 1Z P where 
a relation between microscopic parameters is found, for which the transformation of 
the Hamiltonian into positive semidefinite form holds. 

§3. Exact many-electron ground states on diamond Hubbard chains 




i+r. 



e s = 2 



The first exact many-electron ground 
states for the diamond Hubbard chain 
in an external magnetic field have been 
deduced by Gulacsi, Kampf, and Voll- 
| j>(s =3 /( hardt. 8 ) Here we provide further de- 

tails of the construction process and the 
physical properties of this system. 

The diamond Hubbard chain is 
shown in Fig. 1. This chain has 3 sites 
per unit cell, hence 3 sublattices, de- 
noted by the sublattice index s = 1,2,3, 
and therefore 3 bands. One denotes by 
N, N c , and n = N/(3N C ) the number of 
electrons, cells, and electron density, respectively. The Hamiltonian of the chain 
H = Hq + Hu is given by 



Fig. 1. Diamond Hubbard chain; the hopping 
matrix elements on different bonds are in- 
dicated, as well as the sites on the sublat- 
tices s = 1 and 2 with potential e. 



(J i=l 

+ t±cl +r2 ^c i+ruCT + t||cf +a)<r ci )<r + H.c] n i+rsj(T }, 

3=1,2 

N c 3 
i=l s=l 



(3-1) 
(3-2) 



a) 



b) 




Fig. 2. Schematic view of three particular 
bare band structures of the diamond Hub- 
bard chain for (a) t± = iy = 0, S = it/2, 

(b) i|| > 0, e = -t± + ij 1 /2, 5 = n, and 

(c) t|| =0, et± cos 8 = t± — cos 2 S. 



In the external magnetic field the 
hopping matrix elements acquire the 
Peierls phase factor with 5 = 2ir<P/<Po 
and = hc/e is the flux quantum. 
Here we have chosen the vector poten- 
tial A || a, and the field dependent hop- 
ping amplitudes follow from iy/(B) = 

iy/(0) exp[(2ivr/# ) ff A ■ dl]. For one 
flux quantum per unit cell (triangle) 
8 = 7r. The Zeeman term is not ex- 
plicitly included in H, but all deduced 
ground states have also been analyzed 
in the presence of the Zeeman term. 
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Fourier transforming Hq one finds 
Ho = Sk,cr Es,s'=i -^s,s'(k)4,k,CT<V,k,(T, where the matrix elements M SjS >, which de- 
termine the free-electron dispersion, are given in Ref. 8 ) In the following all energies 
will be given in units of 2t. Examples for bare band structures are schematically- 
shown in Fig. 2. 

3.1. Flat-band ferromagnetism 

We first consider the parameters 



t± 



tn 



0,5 




|, in which case 
the single-electron band structure con- 
sists of three flat bands (see Fig. 2a). 
This corresponds to the "Aharonov- 
Bohm cage" limit analyzed by Vidal et 
al 22 ^ for the two electron problem with 
e = (see also Fig. 3). Based on their 
finding that excited singlet eigenstates 
are localized for U = but delocalized 
for U > 0, Vidal et al. conjectured that 
such delocalized states should emerge 

also at finite electron densities. Below we deduce indeed exact ground states for 
the diamond Hubbard chain at finite densities. 

Following the strategy outlined in §2 one first transforms the kinetic part of the 
Hamiltonian into positive semidefinite form by introducing new canonical fermionic 
operators C v ,\,o m real space (for details see Ref. 8 )); these operators represent lo- 
calized Wannier eigenstates (see Fig. 3). With the energies £2 = £, E2±i = 
(e =F Ve 2 + 4)/2 one obtains 



Fig. 3. Diamond Hubbard chain for t± = t\\ = 
0. The thin line depicts the localized eigen- 
state corresponding to C\ i cr . 



N c 3 
<r,i=l v=\ 



(3-3) 



Since C\ Ka C v ^ C! and Hu are positive semidefinite operators the ground state of H 
for N c > N , U > becomes 



N 



\^m = U d hJ°)> K = E * N - 



(3-4) 



i=l 



For n = 1/3 one obtains a fully saturated ferromagnetic ground state in the 
lowest band (see Fig. 4) which represents an explicit realization of Mielke-Tasaki's 
flat band ferromagnetism. 23 ) The strategy for the proof of the uniqueness of this 
ground state is discussed in Appendix A. For n < 1/3 only the Wannier states with 
a spatial overlap have the same spin, the (highly degenerate) ground state hence 
consists of ferromagnetic clusters of arbitrary spin orientation. Since Hu\$/g(N)) = 
and the kinetic part of H is diagonal in real space, the ground state (3-4) is localized. 
The explicit proof of the uniqueness of (3-4) for n < 1/3 is presented in Appendix 
B. 
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The flat-band ferromagnetism found here realizes early ideas of Gutzwiller 24 ) 
and Kanamori 25 ) from 1963 about the origin of ferromagnetism. They argued that 
the ferromagnetic orientation of the electronic spins is due to the fact that the 
local interaction is thereby completely suppressed thus lowering the energy greatly. 

However, they were not able to com- 
pute the associated loss in kinetic en- 
ergy of this highly correlated state. In 
the present case of flat-band ferromag- 
netism the latter effect poses no problem 
since the kinetic energy is zero anyway. 

We note that at U > the lowest 
band is flat only for S = vr/2, but is dis- 
persive for 5 = it, or 5 = 0. In the latter 
situation the system is most probably 
conducting for incommensurate fillings. 
Hence one encounters here a magnetic 
field induced metal-insulator transition. 

3.2. Correlated half-metal 

The previous results suggest that 
itinerant states are easier to realize at 
5 it/2. To analyze this problem we in- 
vestigate below the special case tu > 0, 
e = — 1± + t J ^ 1 /2, and 8 = it; the single- 
electron band structure has a dispersive 
lowest band and two upper flat bands (see Fig. 2b). 

3.2.1. Transformation of H into positive semidefinite form 

For the rewriting of Hq one defines first the non-canonical fermionic operators 

A\,o = dicker + a2Ci+r 2 ,cr + a?,C\.+&,CT + 0,4^+^,(7 , (3-5) 

which fulfill (i i>(T ) 2 = and {i i>tr , ij j / fly. Based on (3-5) one obtains 
A \,a A h° = (aWci+r^^ + a* 3 a 2 cl +ai(T c i+r2ta + a* 4 a 3 c\ +rua c i+ata + a\a A c{ a c i+r ^ a 

+ a2«4Ci+r 2 ,<7 a i+r 1 , ( 7 + a3 a lCj[+a,<A<r + H - C ') 

+ \a\\ 2 n it(T + |a 2 | 2 n i+r2)(T + |a 3 | 2 n i+aj(T + |a 4 | 2 n i+rij(T . (3-6) 

As discussed in Sec. 2.2. the coefficients oi, a 2 , 03, 04 can be determined in different 
ways, namely by either requiring -Ei, ff 4,A" = ^ or + A,*^ = ^o- 
Here we choose the former condition which will enable us to construct exact ground 
states for N > 4N C (see Sec. 3.2.2.). With the requirement - £\ A\ a A^ a = H 
one arrives at the matching conditions 

a 2 a\ = a%a 2 = a\a 3 = a^a^ = —te 1 ?, a 2 a 4 = —t±, a^a\ = —t\\, 

I a-i | 2 + |o-3| 2 = e + |a 2 | 2 = e + |a 4 | 2 , (3-7) 




k 

Fig. 4. Ferromagnetic flat-band solution on 
the diamond Hubbard chain obtained for 
t± — £|| = 0; the ground-state result does 
not provide information about higher lying 
states as indicated by the question mark. 
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from which 

A,a = y/t\\[Ci,a- ~ Ci+a,<r ~ 2t_ L e % 2(c i+rijCT - C\+t 2 ,it)\ (3-8) 

follows. Based on (3-6 -3-8), Hq is rewritten as 

4 

i,cr i,cr m=l 

For the interaction part Hy, one uses 

3AT C 

H v = C/^n iiT n u = UP + UN- UN C , P = ^P;. (3-10) 

i i 

iV is the particle number operator, and Pi = (n^f — — 1) is a positive semidef- 

inite operator, which gives zero for at least one electron on site i and 1 for an unoc- 
cupied site i. Consequently, H has been transformed into the positive semidefinite 
form 

H = K A {c + Up + 4 7 > ( 3 - n ) 

i,er 

where E 1 / = (e + U + t ± )N - N C (3U + 4t± + l/t ± ). 
3.2.2. Construction of the ground state 

As discussed in Sec. 2.2, the ground state \^ g ) for U > is deduced such to 
satisfy -Ajj^g) =0 and P\& g ) = 0, from where indeed 

H\* g ) = E g \V g ) (3-12) 
follows. Since (i| ff ) 2 = one finds 

1^(4^)) oc II K-J{A*) ■ (3-13) 

i 

In order to satisfy P\& g ) = as well, which requires one electron on each site, besides 
the contribution from (3-13) one has to introduce in the ground state wavefunction 
also the operator 

P t = Ul4 + r Sia A + r Si ,J> (3-14) 
i 

which creates two electrons with fixed spin a on arbitrary sites of each unit cell. 
Hence, the unnormalized ground state becomes 

|Sf (42V C )> = ciHAl^AljF^O), (3-15) 

i 

where c is a normalization constant. The contribution in (3-13) creates one a elec- 
tron, and an electron with spin —a in each unit cell, while P\ creates two electrons 
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with spin a on arbitrary sites of each unit cell. Consequently the ground state in 
(3-15) corresponds to N = 4N C , n = 4/3, n a = N a /3N C = 1, n_ CT = 1/3. The proof 
of the uniqueness of (3-15) is presented in detail in Appendix A. In conclusion, (3-15) 
is the unique ground state. 

Concerning the physical properties of the ground state one observes that in 
(3T5) each lattice site is occupied by one cr-electron; hence, the cr-electrons are 
localized. The —a electrons are spatially extended, but localized for N c — ► oo. This 
is inferred from the ground-state expectation value of the hopping term r r - a = 
(c. _ CT cj +r _ ct + H.c), which in the limit N c — ► oo, for r/a = m, is calculated as 

The one-particle localization length £_ CT in (3-16) increases almost linearly with 
l/t±- 8) 

3.2.3. Solution for n > 4/3 

For electron densities n > 4/3 there are N > 4N C electrons. No cr-electrons but 
AN — cr-electrons can be added to the system such that n a = 1, n_ CT = 1/3+ AN/N C . 
The ground state then becomes 

AN 

|Sf (4JV C + AN)) = [J 4«, ka ,-J^(^c)), (3-17) 

a=l 

where n a can be arbitrarily chosen from s = 1,2,3. Eq. (3-17) contains plane 
wave-type contributions from the — cr-electrons. Consequently, the additional AN 
—a electrons are itinerant. Indeed, this follows explicit from the expectation value 
of the hopping term r r ^ a . For |r/a| S> 1, AN = 1, s = 1, one finds r r a / r^_ a = 

1 — [2 + (1 — cosa/c)/t J -1 . Here i^-er ^ s * ne pl ane wave result and k is the norm 
of the momentum of the electron added above n = 4/3. We note that the momenta 
k Q in (3-17) are arbitrary. Hence it is not possible to define a Fermi surface, and the 
ground state therefore corresponds to a non-Fermi liquid. 

In the density range 4/3 < n < 5/3 the ground state has 3A^ C immobile a 
electrons and N c electrons with spin —a confined to their localized Wannier functions. 
For this reason the ground state has a high geometric degeneracy; its uniqueness is 
discussed in Appendix B. 

Only the AN conducting —a electrons are itinerant, leading to metallic behavior 
with a low carrier density and a low spin polarization. Since the conduction through 
this correlated half-metal involves only electrons of one spin species, such a system 
may serve as a spin-valve device. Since U > is required for this property one 
encounters here a correlation induced localization-delocalization transition to a half- 
metal. The magnetization behaves as M oc (1 — AN/N C ) — > for AN — ► N c . 
in a finite magnetic field. The transition can also be induced by tuning the local 
potential e. At n = 5/3 the ground state becomes nonmagnetic and localized in 
the thermodynamic limit, while at n > 5/3 one finds a nonmagnetic but conducting 
ground state. 
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3.3. Exact ground states for general magnetic flux 

It is also possible to construct exact ground states for general magnetic flux 
values including also the B = case, for example at 5 € (— 7r/2,7r/2) with tn = 0, 
t± < 0, b = — cos 6/t±, e = b — b^ 1 , and n > 5/3; 8 ) the single-electron band structure 
has two dispersive bands and an upper flat band (see Fig. 2c). For these hopping 
and potential parameters one obtains a localized nonmagnetic ground state for 5 = 
over an extended range of densities n > 5/3. By contrast, at finite magnetic fields 
the ground state is a nonsaturated ferromagnet. This localized state at n = 5/3 is 
insulating, but gapless and itinerant for n > 5/3. For the latter densities the majority 
spin cj-electrons are immobile and only the —a electrons are itinerant. Therefore, by 
varying the external magnetic field and the sublattice potential one can tune from 
a localized, nonmagnetic ground state to a nonsaturated ferromagnet in the density 
range n > 5/3; the ferromagnetic state is insulating at n = 5/3, but conducting for 
n > 5/3. 

We finally note the properties derived above do not depend on the Zeeman 
coupling. Only the solution in Sec. 3.1 for n < 1/3 is altered by the Zeeman term, 
whose presence leads to a fully spin aligned ferromagnetic solution. 

§4. Exact many-electron ground states on triangle Hubbard chains 

The procedure for the construction of exact ground states on the diamond Hub- 
bard chain is now applied to a Hubbard chain composed of corner-sharing triangles 
(see Fig. 5). Triangles are well-known to frustrate non- ferromagnetic magnetic or- 
der. Hence interacting electrons on lattices with triangles as structural sub-units are 
expected to favor ferromagnetic ground states. 

The triangle Hubbard chain in Fig. 
5 has two sites per cell which leads 
to two energy bands. The number 
of electrons, cells, and electron den- 
sity are here denoted by N, N c and 

n = N/(2N C ), respectively. This model J+f jj 

was previously investigated in detail by 1 
Penc et al. 2 ^ Fazekas et al. 27 ^ Miiller- Fig. 5. Triangle Hubbard chain 

Hartmann 28 ) and Derzhko et al. 29 ^ with 
a variety of analytical and numerical techniques. 

The Hamiltonian of the triangle chain H tri = H^" 1 + Hy l has the form 

N c 




= E E [ (4* W + *'4^+ r >- + 4 +a ,A+r,a) + H.C.) 
a i=l 

Hfri = Y^[U n iA h id + C/in i+r)T n i+r J. (4-1) 



The model parameters are explained in Fig. 5. In the following we will construct 
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exact ground states in the subspace of parameters defined by the relation (t') 2 /t = 
£1 — £o + 2£, E\ — eo > —2t, t > 0, for which the lowest bare band is flat. 

4.1. Ground states for the triangle Hubbard chain 
4.1.1. Ground state for C/ , f/i > 



a) 




b) 




Fig. 6. Ground-state configuration on the triangle Hubbard chain for (a) Uo,Ui > 0, and (b) 
Uo > 0, Ui — 0. The ground-state results do not provide information about higher lying states 
as indicated by the question marks. 



Employing the operators 

A,a = Vt[ h,a + Ci+a,a + (t' /t)c i+T:a ] (4-2) 

the Hamiltonian (4-1) can be transformed into positive semidefinite form: 

fitri = A\ a A i>l7 + KN + H%\ (4-3) 

i,cr 

Here N is the particle number operator and K = eq — 2t. To construct the ground 
state one has to find operators B[ a which satisfy {A- hcr ,Bj, a ,} = for all values of 
the indices. The ground state then has the form 

N 

\V tri (N)) = l[Bl i \0}, (4-4) 

i 

with b\ = [c|_a+r,<7 + d l+r,a ~ C*'/*)^], N < N c . For n < 1/2 there exist ferro- 
magnetic clusters, and only operators b\ which touch have the same spin index. 
For n = 1/2 all clusters touch such that a fully saturated ferromagnetic state is 
obtained as shown in Fig. 6a for the Ui > case. The ground state corresponds to 
a particular realization of flat-band ferromagnetism, 23 ) for which Derzhko et al. 29 ^ 
calculated the thermodynamic properties. 
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4.1.2. Ground state for U > 0, Ui = 

The ground state has the same form as (4-4), but now it is non-magnetic even 
for n = 1/2 (see Fig. 6b). This is due to the fact that, although the connectivity 
condition is fulfilled (i.e., Bj a and Bj, a(j overlap at site i + r) there is no interaction 
at the connectivity point (Ui = 0). The fulfillment of the connectivity condition for 
the Wannier functions is therefore not sufficient for the ground state to be ferromag- 
netic. In addition a finite Hubbard repulsion is needed at the connectivity points to 
guarantee a correlation between the electrons. 

4.1.3. Ground state for U = 0, f7i > 

The ground state has again the form (4-4), but now it is a fully saturated ferro- 
magnetic state at n = 1/2 even for Uq = 0, since the Hubbard interaction is non-zero 
at the connectivity points (U\ > 0). We thus see that for flat-band ferromagnetism 
to occur the Hubbard interaction is needed only at the connectivity points. The 
ferromagnetic states discussed in Sec. 4.1.1 and Sec. 4.1.3. are localized in the 
thermodynamic limit. 

4.2. Mapping to the one- dimensional periodic Anderson model (PAM) 
By the following change of notation in (4-1) 

ha -> Ci+r, CT -> fi,<7, t' — ► V, E X -> Ef, Ul -> U, (4-5) 

and setting eq = Uq = the triangle Hubbard chain transforms into a one-dimensional 
PAM with zero on-site hybridization and nearest-neighbor hybridization V: 

N c 
a i=l 

+ V(4J i+T , a + d\ +a J i+r , a + H.c.)] "f., (4-6) 

i 

A lowest bare flat band is obtained for V 2 /t = Ef + 2t, Ef > -2t, t > 0. 

4.2.1. Exact ground state for the one-dimensional PAM 

In the parameter regime V 2 /t = Ef + 2t, Ef > — 2t, t > the exact ground 
state for the one-dimensional PAM is obtained from (4-4) as 

N c 

w pam {n = Nc)) = H { fl a+r a + ^ _ — 4J|0>. (4-7) 

i=l 

4.2.2. Itinerant ferromagnetism of finite chains 

In analogy with (3T6) the localization length £ of the (i-electrons can be cal- 
culated from the expression for the expectation value of the hopping term r r . Its 
dependence on the nearest-neighbor hybridization strength V is given by 
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as shown in Fig. 7. Here a is the distance between sites at the base of the triangle 
chain (see Fig. 5). Although flat-band ferromagnetism in the thermodynamic limit 
usually corresponds to localized electrons, the localization length can be so large 
that it exceeds the length of a finite sample. In realistic situations the exact ground 
state of the finite chain therefore shows itinerant ferromagnetism. 

4.3. Application to f -electron ferromagnetism in CeRh^B2 

As pointed out by Kono and Kuramoto 21 ) the one-dimensional PAM, (4-6), 
may be employed to model the 4/-electron compound CeRh 3 B 2 . This material is 
ferromagnetic below T c = 120 K. Its properties are interesting for several reasons: 
(i) the ferromagnetic state cannot be explained by a RKKY interaction, (ii) the 
4/ moment has the remarkably small value of 0.45 hb compared to the free Ce 3+ 
moment of 2.14 /ie, and (iii) the Curie temperature T c has the highest value among 
all known Ce compounds with non-magnetic elements. 

The lattice structure of the class of 
rare-earth (RE) compounds RERI13B2 £ 
was determined by Yamada et al. 30 ^ 
A cut through the (a,c)-plane of the 
Cerium system (RE = Ce) is shown in 
Fig. 8. The remarkable difference be- 
tween the distances of the Ce atoms 
along the a- and the c-direction, respec- 
tively, makes this material a quasi one- 
dimensional two-band system. 31 ) For 4 ^ ^ ^ ^ ^ 

this reason Kono and Kuramoto 21 ) mod- j 

eled CeRh 3 B 2 by the triangle PAM 

, . , • -n- n mi Fie. 7. Localization length £ of the d-electrons 

chain shown m .big. 9. Ihe parame- ° il ,. . , „ N 

of the one-dimensional PAM, (4-6), as a 

ters entering the model are t = 0.34 eV, function of the nearest-neighbor hybridiza- 

V = 0.24 eV, E f = -0.714 eV, U = 7 tion strength V. 

eV, n = 0.55, as obtained from band- 
structure calculations and X-ray absorption spectroscopy (we note that the density 
n = 0.55 corresponds to twice the value, i.e., n = 1.1, in the notation of Ref. 21 )). 

4.3.1. Variational results for the one-dimensional PAM 

Kono and Kuramoto 21 ) computed the band structure of CeRh3E>2 by construct- 
ing a Gutzwiller variational wave function and using variational Monte Carlo to 
determine the ground state. They found the lowest band of CeRli3B2 to be almost 
perfectly flat and fully polarized, and the second band dispersive and almost un- 
polarized. It should be stressed that the shape of the lowest quasiparticle band of 
interacting electrons depends on U. 

4.3.2. Exact ground state of the one-dimensional PAM 

We can use the exact ground state of the one-dimensional PAM, (4-7), to make 
direct contact with the /-electron ferromagnetism found in CeRh3B2. Namely, this 
ground state represents a saturated flat-band ferromagnetic state of the /-electrons 
for electron density n = 1/2 and arbitrary U > 0, provided the condition V 2 /t = 
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Fig. 8. Cut through the (a,c)-plane of the lat- 
tice structure of CeRhaB2 based on the pa- 
rameters determined in Ref. 30 ' 
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Fig. 9. One-dimensional PAM used to model 
CeRh 3 B 2 . 21 > 



Ej + 2t, Ef > —2t, t > holds. The condition is indeed fulfilled for the parameter 
values t = 0.34 eV, V = 0.23 eV, E f = -0.52 eV, which are very close to the 
experimentally determined values. In this case the lowest band is perfectly flat for 
all U > 0. Apparently the ferromagnetism of the /-electrons in CCRI13B2 is due to 
the presence of an (almost) flat lowest band. 

4.3.3. Magnetic moments 

The experimentally observed magnetic moment of the /-electrons is rrif = 0.45 
/i#. 32 ) The variational result 21 ) yields rtif = 0.94 [1b, while the exact ground state 
discussed above leads to a value of m.f = 0.68 hb which is even closer to the experi- 
mental value. Both theoretical values are somewhat higher than the experimental / 
moment. This may be due to the fact that the theoretical approaches do not include 
the spin-orbit interaction which causes an antiferromagnetic coupling between the 
d- and the /-electrons. 



§5. Conclusions 



Despite its simplicity, the diamond Hubbard chain displays remarkably complex 
physical properties with insulating, conducting, and fully or partially spin polarized 
states. The constructed exact many-electron ground states revealed flat-band fer- 
romagnetism, correlated half-metallic behavior, and metal-insulator transitions in 
the diamond Hubbard chain. The virtue of tuning fundamentally different ground 
states through external magnetic fields, site selective potentials, or the electron den- 
sity, points to new prospects for designing electronic devices, which allow to switch 
between different ground states, including states that operate as a spin valve. We ap- 
plied the general strategy for the construction of exact ground states also to triangle 
Hubbard chains. The mapping of the triangle Hubbard chain to the one-dimensional 
periodic Anderson model within a restricted parameter set provided a theory for 
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/-electron ferromagnetism in CeRli3B2, which is based on the physical properties of 
an exact ground state. 

Further applications of the construction of exact ground states by casting the 
Hamiltonian into positive semidefinite form are possible for model systems with 
selected geometries in artificially designed structures. In particular, we constructed 
exact ground states for pentagon Hubbard chains. These results will be presented 
elsewhere. 
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Appendix A 

Uniqueness of ground states with N = pN c where p is an integer 

Here we prove the uniqueness of those exact ground states constructed in this 
paper for which the electron number N is an integer multiple of the number of unit 
cells N c . The proof proceeds as discussed in Sec. 2.3 and is exemplified by proving 
the uniqueness of the ground state l&g 1 (4N C )) , (3T5), of the diamond Hubbard chain. 
In this case N = 4iV c , i.e., p = 4, and the electron density has the value n = 4/3. The 
Hamiltonian, the unnormalized ground state, and the A\ a operators entering therein 
are given by (3-11), (3T5), and (3-8), respectively. To prove the uniqueness one first 
has to determine the kernel of H' = H - E l g l = Y, c Ya=\ ^ A \a + u ?- The form 
of the Hamiltonian H' implies that its kernel is simply given by the intersection of 
the kernels of the individual terms, i.e., ker(_£T) = ker(^ CT Y^a=i A,(r^L) fl ker(P), 
where P is an operator, which assumes its minimum eigenvalue 0, if there is at least 
one electron at each site. 

To prove the uniqueness of the ground state l&g 1 (4N C )) , (3-15), it is now neces- 
sary to demonstrate that it spans ker(H'). According to the discussion in Sec. 2.3 
the proof consists of two parts. First one has to establish that 

(i) l&g 1 (4:N C )} is an element of the intersection of the kernels of the positive 

semidefinite operators entering in H' . 
Then one must show that 

(ii) all states from the kernel of H' can be written in the form of the ground 
state |^ J (47V C )). 

The first part of the proof, (i), is covered by the following Theorems A-l to A-3. 
The second part of the proof, (ii), then follows immediately (see below). 

To prove (i) we note that the first term in H' is a sum of positive semidefinite 
operators Ai :a A\ . The ground state of H' must therefore be an element of each 
kernel of all the operator Ai jCr Aj a , for all i, a. Thus we first determine the kernel of 
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(Sec. A.I.), then the kernel of the sum E CT Eil'i A^i l(J ( Sec - A - 2 -) 5 and 
finally the kernel of the full Hamiltonian H' (Sec. A. 3.). 

A.l. The kernel of Ai^Al^ 

Theorem A-l: The kernel of Ai )(T Aj a is spanned by states of the form 

|^) = 4>t| ), (A.l) 

where is an arbitrary operator. 
Proof of Theorem A-l 

a) From (if j 2 = it follows that A ij<T A\ \&) = 0. This implies that \&) G 
kev(A iia A[ a ). 

b) In order to show that all elements of ker(Ai iCr Aj a ) can be written in the form 
\&) = A\^W^\0) we consider a state \<P) = Y^\0) G kev(A ija A\ a ) and show that 
has the form |#) = A^W^O). 

Making use of {A i>er , A\^ = 8t 2 ± t\\ + 2t\\ =: a? ; > and A^ a A\^\<P) = we can 
write 

I*) = 5*K» = ^(Aai + AU^IO) = Al^A^M. (A-2) 

"i,i "i,i 

With the identification := -^-A\ a Y^ we have therefore verified that indeed 

a i,i ' 

W) = A{ a W^) = W). (A-3) 
This completes the proof of Theorem A-l. 
A.2. The kernel of £^ M^\, 

The kernel of Ylu^i ^i^^L has the form 

N c N c 

me e KKr) = n n m^mu- ( a - 4 ) 

i=l <r=T,ii=l 
Theorem A-2: ker(^ CT Yli^i ^i,crA\ a ) is spanned by vectors of the form 

N c 

cr i=l 

where the operator W' is arbitrary. 

Theorem A-2 follows directly from the linear independence of the operators 
Aj CT (that is, Aj ct / Sj^i a j,o-^j a holds for arbitrary coefficients a- hC due to the fact 

that the operators cj, „ j cj, „ enter only in A? and not in At , for j 7^ i. 

i+ri,<r' i+r2,cr J i,cr j,<t' J ' 
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A. 3. The kernel of H' 

The kernel of H' is given by ker(H') = ker(^ CT Ylu^i a ) fl ker(P), where 

P is an operator, which assumes its minimum eigenvalue 0, if there is at least one 
electron at each site. Hence, ker(H') is spanned by vectors of the form in Eq. (A-l), 
where now the operator is no longer arbitrary, but must have a specific form to 
satisfy the abovementioned requirement imposed by ker(P). 

Theorem A-3: For N = 4N C the operator is given by = F<J, where F„ 
creates two electrons with fixed spin a in each unit cell. 
Proof of Theorem A-3. 

The proof proceeds in two main steps: we show that = F„ is 1) a possible 
choice, and 2) the unique choice. 

1) Let us denote by Ti the unit cell (= triangle), which contains the sites i, i + 
ri, i + T2- The state F,J|0) has one empty site in each unit cell Ti. We now analyze 
how these empty sites are filled by the product Yls^i ^L- We first observe that A^ 
acts on all sites in Tj, which implies that the empty site in Ti can indeed be filled 
by applying A\ a . There are N c empty sites and N c operators A\ a \ thus, the product 

n£i^iCT inserts N c electrons with spin a into the system. Since (cj CT ) 2 = 0, the 
state njI c i^L£Jl°> has 

one cr-electron on all 3N C sites of the lattice. Hence, Fj is 
a possible choice for in constructing the states contained in ker(i7'). 

2) We now flip the spin a of one electron in F„ into —a, and denote the new 

operator W' by P' a . In this case, a single unit cell, say Ti o , will have one empty site, 
one electron with spin a on another site, and one electron with spin — a on the third 
site. Then the operator A] _ can create one electron with spin a on the site, which 
already contains the electron with spin —a, thereby creating a double occupancy. 
One empty site will remain in Ti o . For i ^ i D , the operators Aj can introduce their 
electron into the empty sites of Ti. Hence one empty site can remain in the system. 

The state |x) = Ili^i a F'a\ty therefore has at least one term with an empty site 
and contains one —a electron. Acting now with fliti M-a on lx) adds N c electrons 
with spin —a on (3N C — 1) sites. The new state \<P) = Yii^i A\ _Jx) wu 1 therefore 
have at least one term with an empty site, and is thus not an element of ker(P) and 
ker(iJ). This shows that = PI is the unique choice. This completes the proof of 
Theorem A-3. 

Theorems A-l, A-2 and A-3 prove that |l^J / (4A c )), (3-15), is indeed an element 
of the intersection of the kernels of the positive semidefinite operators entering in 
H' . This completes part (i) of the proof of the uniqueness discussed in Sec. 2.3. 

Since all vectors from ker(H') can be written in the form |^ 7 (4A C )), (3-15), and 
since this state represents the constructed ground state 1^}, it follows immediately 
that all vectors in ker(fr) can be written in the form (3-15). Therefore the ground 
state (3-15) is the unique ground state of the Hamiltonian H, (3-11). 

The proofs of the uniqueness of the ground state (3-4) at density n = 1/3, 
corresponding to p = 1, and of solution III in Ref. 8 -* at density n = 5/3 (p = 5) 
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proceed similarly. 

Appendix B 

Uniqueness of ground-state solutions with N = pN c for non-integer p 

B.l. Geometric degeneracies 

For non-integer p as in the case of (3-4) for n < 1/3, and (3-17) for n > 4/3 (as 
well as in the case of solution III for n > 5/3 of Ref.p]), the ground states have a high 
geometric degeneracy. In these cases there exist linearly independent eigenstates \^g) 
with the same ground-state energy. They may be denoted by a degeneracy index 
7 = 1,2, ...,Md, where is the degree of degeneracy. The ground state then 
becomes \& g ) = ^7^1 a j Wg)i where apart from the normalization condition, the 
coefficients a 7 are arbitrary. For these ground states the proof of uniqueness proceeds 
analogously to that described in Appendix A. Namely, one has to prove that \^P g ) 
spans kev(H'), but now, given the specific structure of \& g ), it has to be verified that 

1) all linearly independent states \^g) are elements of the intersection of the kernels 

L 

of the positive semidefinite operators P n , i.e., \\Pg) € ker(H') = f) ker(P„), and 

n=l 

L 

2) all states 6 f] ker(P n ) can be written in the form of the constructed ground 

71=1 

state, e.g. as a linear combination of the vectors \&g). 

As an example we now prove the uniqueness of the ground state \&g(N)}, (3-4), 
for n < 1/3. The proof is presented for fixed N, namely N = N c — 1. For other N 
values the proof proceeds similarly. The first part of the proof, (i), (see Sec. 2.3) is 
covered by the following Theorems B-l and B-2. The second part of the proof, (ii), 
then follows immediately (see below). 

B.2. Uniqueness of the ground state (3-4) for N = N c — 1 particles 

The transformed Hamiltonian H = Hq + Hy is determined by Hq and Hu as 
given by (3-3)and (3-2), respectively. Furthermore, one has iV = J2\ a Y%=\ i (fiv,\,ai 
where C v ^ a are genuine canonical Fermi operators, and E% < E\,Ei- Since iV is a 
constant of motion, the Hamiltonian can be written as H = Y2i CT [(-^i — ^3)^1 i (J C'i,i,<T+ 
(E 2 - E z )Cl Xa C 2X(T ] + H V + E\, where E\ = E 3 N and {E l - E 3 ), (E 2 - E 3 ) > 0. 
Consequently, H' := H — Eg is a positive semidefinite operator: H' = P\ + P 2 , where 

A = Ei, CT [(^i " E 3 )C\. ;a C 1X(7 + (E 2 - E 3 )Cl. !a C 2X(J ], and P 2 = H v . The ground 
state is given by (3-4). 

B.2.1. The kernel of A 

Theorem B-l: The kernel of P\ is spanned by states of the form 



N c -1 

i^(o},w)>= n ^i°>. 

i=l 



(B-l) 
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where the product is taken over arbitrary N c — 1 sites, and dj is arbitrary. Note that 
|^({i}, {cJi})) depends on the set of sites {i}, and the set of spin orientations {ci}. 
Proof of Theorem B-l. 

The proof proceeds in two steps: a) one first shows that Eq. (B-l) is a possible 
choice, and b) then verifies that it is the unique choice. 

a) The initial 6N C canonical Fermi operators Ci +rsi<J , where s = 1,2,3 is the 
sublattice index, have been transformed by a linear transformation into the 6iV c 
opertors C v i ai , which satisfy canonical fermionic anti-commutation rules. Since the 
set of the new operators {C^i,^} substitutes the set {ci+ rs)Cr }, all operators in the 
system acting on states defined by arbitrary (i, a\) can be expressed in terms of the 
operators C U;i>ai . 

Given these conditions, an operator X^ which satisfies 

{C 1Xav X^ = {C 2;hai ,X^ = (B-2) 

for all indices i, o\, must have the form X^ = T\- T\ ct ■ . Indeed, the anti- 

commutation relations (B-2) hold only for this type of X^ operator since an arbi- 
trary Xt operator can be expressed in terms of C u j >a . operators, and only C^. ja , 

anticommutes with all Ci,i,o-, and C2,i,a r 

b) If X^ contains at least one C^- operator with v = 1,2, the relation Eq. 
(B-2) no longer holds, hence Eq. (B-l) is the unique choice. This proves Theorem 
B.l. 

We furthermore observe that due to the fermionic operator nature of C3 i (j,, the 
deduced elements of ker(Pi) from Eq. (B-l) are not only linearly independent but 
also orthogonal: 

<#({i}, {o-ijwm, W)) = <%,{i'}<^},w}- ( B -3) 

B.2.2. The kernel of H' 

Since ker(i?')=ker(Pi) n ker(P2) 5 the states in ker(P') have the form (B-l), but 
the specific set of indices ({i}, {cii}) has to be chosen such that these states are also 
contained in ker(p2)- This condition requires the proof of the following theorem. 

Theorem B-2: The kernel of H' is spanned by states of the form 

iv c -i 

|ff(j)> = l*({I})> = II ^|0>, (B-4) 

i=l 

where a is fixed, and the product extends over N c — 1 arbitrary sites. 

Since the product in Eq. (B-4) is taken over iV c — 1 sites of the sublattice s = 3, 
the vector |^({i})) can be denoted by the specific site j excluded from the product, 
i.e. (^({i})) = |^(j)). This site j plays the role of the degeneracy index 7 in Sec. B.l. 

Proof of Theorem B-2. 
The theorem is proved in three steps: First one has to show that all a\ spin indices 
in (B-4) must have the same value, secondly it has to be verified that the deduced 
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states in ker(H') are linearly independent, and thirdly one needs to prove that the 
set of the deduced states is complete. 

a) All states in ker(P2) require the absence of any double occupancy. Therefore, 
all indices a\ must be fixed to the same value a. The double occupancy cannot be 
avoided otherwise, since in N c unit cells iV c — 1 operators Cg j share a common site. 

b) due to the relation (B-3) the vectors \&Qi)) are linearly independent, since 
they are orthogonal = 5j y. 

c) Obviously, the degree of degeneracy of a macroscopic many-body state is given 
by the number of equivalent configurations. Due to part a), apart from the trivial 
degeneracy with respect to the total spin direction, on finds = N c , because there 
are N c possibilities for placing N c — 1 identical objects on N c sites. Consequently, 
ker(H') is an = N c dimensional Hilbert subspace. The index j in the states of the 
form (B-4) refers to N c possible different unit cells. Furthermore, the vectors |^(j)) 
are linearly independent, hence the set {|^(j))} forms a (orthonormalized) basis of 
ker(H'). Therefore, {|^(j))} spans ker(H'), and the ground state has the general 
form \\P g ) = J^j^i with arbitrary coefficients ay This proves Theorem B-2. 

Since {|^(j))} is a basis of ker(H'), all states \<f>) 6 ker(H') can be written in the 
form |<Pg) of part c). Therefore (3-4) is the unique ground state of the Hamiltonian 
H consisting of the two terms (3-3) and (3-2). 

For N = N c the set {j} is an empty set, and the upper limit set to N c in Eq. 
(B-4) provides the unique ground state at density n = 1/3 (see (3-4) for N = N c ). 
A similar strategy for N < N c — 1 allows one to prove the uniqueness of the ground- 
state solutions for all electron densities n < 1/3 in (3-4). Similar proofs are expected 
to hold in the case of the ground state |^ 7 (4iV c + AN)), (3-17), for n > 4/3 (p > 4) 
and solution III for n > 5/3 (p > 5) constructed in Ref. 8 ^ 
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